#### EXERCCIOS 5 - ANA CAROLINA PAVAN #####

#### Editando alguns parmetros grficos ####
esaligna <- read.table("esaligna.csv", header=TRUE, sep=",", dec=".")
str(esaligna)
#1
plot(dap~ht, data=esaligna, xlab="DAP(cm)", ylab="Altura(m)")
#2
plot(dap~ht, data=esaligna, tcl=0.3, xlab="DAP(cm)", ylab="Altura(m)")
#3
plot(dap~ht, data=esaligna, tcl=0.3, bty="l",xlab="DAP(cm)", ylab="Altura(m)")
#4
plot(dap~ht, data=esaligna, tcl=0.3, bty="l", xlab="DAP(cm)", ylab="Altura(m)", main="DAP x Ht")
#5
plot(dap~ht, data=esaligna, tcl=0.3, bty="l", cex=2, xlab="DAP(cm)", ylab="Altura(m)", main="DAP x Ht")

#### Dois grficos juntos ####
esaligna<-read.table("esaligna.csv", header=TRUE, sep=",", dec=".")
str(esaligna)
esaligna1 = esaligna[,3:4]
str(esaligna1)
as.factor(esaligna1$talhao)
media.dap <- aggregate(esaligna[4],by=esaligna[3],FUN=mean)
desvio.dap <- aggregate(esaligna[4],by=esaligna[3],FUN=sd)
media.dap$desvio<-desvio.dap$dap
media.dap
str(media.dap)
esaligna_medias<- as.matrix(media.dap)
class(esaligna_medias)
par(mfrow=c(1,2))
boxplot(dap~talhao, data=esaligna1,ylab="DAP(cm)", xlab="Talho")
plot(dap~talhao, data=media.dap, xlim=c(0,6.5), ylim=c(6,21), xlab="Talho", ylab="Mdia+desvio-DAP")
points (x=1,y=media.dap[1,2], pch=16)
points (x=2,y=media.dap[2,2], pch=16)
points (x=3,y=media.dap[3,2], pch=16)
points (x=4,y=media.dap[4,2], pch=16)
points (x=5,y=media.dap[5,2], pch=16)
points (x=6,y=media.dap[6,2], pch=16)
rect(xleft=1,ybottom=(media.dap[1,2]-media.dap[1,3]), xright=1, ytop=(media.dap[1,2]+media.dap[1,3]))
rect(xleft=2,ybottom=(media.dap[2,2]-media.dap[2,3]), xright=2, ytop=(media.dap[2,2]+media.dap[2,3]))
rect(xleft=3,ybottom=(media.dap[3,2]-media.dap[3,3]), xright=3, ytop=(media.dap[3,2]+media.dap[3,3]))
rect(xleft=4,ybottom=(media.dap[4,2]-media.dap[4,3]), xright=4, ytop=(media.dap[4,2]+media.dap[4,3]))
rect(xleft=5,ybottom=(media.dap[5,2]-media.dap[5,3]), xright=5, ytop=(media.dap[5,2]+media.dap[5,3]))
rect(xleft=6,ybottom=(media.dap[6,2]-media.dap[6,3]), xright=6, ytop=(media.dap[6,2]+media.dap[6,3]))

#### Adivinhando o cdigo ####
tabela<- read.table("exercicio3.csv", header=TRUE, sep=",")
tabela
par(mfrow=c(1,2),bty="l", tcl="0.3")
plot(y1~x1,data=tabela, xlim=c(0.5,2.2),ylim=c(0.0,3.0), xlab="Log(Patch size)(ha)",ylab="Euclidean distances",pch=17)
sort(predict(lm(tabela$y1~tabela$x1)))
segments(x0=0.5,y0=2.2,x1=2.2,y1=1.17,lwd=2)
mtext("a", side=3, line=1,cex=2, adj=1,family="serif")
boxplot(x2~y2,data=tabela,ylim=c(0.0,3.0),names=c("Small","Medium\nEdge","Medium\nInterior","Large\nEdge","Large\nInterior","Control"),cex.axis=0.7,outline=FALSE)
par(new=TRUE)
text(x=c(1:5),y=3.0,c("*","*","**","*","***"),cex=1.5)
mtext("b", side=3, line=1,cex=2, adj=1,family="serif")

#### FIM! ####












